Bulk, surface and corner free energy series for the chromatic polynomial on the square and triangular lattices

We present an efficient algorithm for computing the partition function of the q-colouring problem (chromatic polynomial) on regular two-dimensional lattice strips. Our construction involves writing the transfer matrix as a product of sparse matrices, each of dimension ~ 3^m, where m is the number of lattice spacings across the strip. As a specific application, we obtain the large-q series of the bulk, surface and corner free energies of the chromatic polynomial. This extends the existing series for the square lattice by 32 terms, to order q^{-79}. On the triangular lattice, we verify Baxter's analytical expression for the bulk free energy (to order q^{-40}), and we are able to conjecture exact product formulae for the surface and corner free energies.Comment: 17 pages. Version 2: added 4 further term to the serie

Dynamic Critical Behavior of the Swendsen-Wang Algorithm: The Two-Dimensional 3-State Potts Model Revisited

We have performed a high-precision Monte Carlo study of the dynamic critical behavior of the Swendsen-Wang algorithm for the two-dimensional 3-state Potts model. We find that the Li-Sokal bound ($\tau_{int,E} \geq const \times C_H$) is almost but not quite sharp. The ratio $\tau_{int,E} / C_H$ seems to diverge either as a small power ($\approx 0.08$) or as a logarithm.Comment: 35 pages including 3 figures. Self-unpacking file containing the LaTeX file, the needed macros (epsf.sty, indent.sty, subeqnarray.sty, and eqsection.sty) and the 3 Postscript figures. Revised version fixes a normalization error in \xi (with many thanks to Wolfhard Janke for finding the error!). To be published in J. Stat. Phys. 87, no. 1/2 (April 1997

On the non-ergodicity of the Swendsen-Wang-Kotecky algorithm on the kagome lattice

We study the properties of the Wang-Swendsen-Kotecky cluster Monte Carlo algorithm for simulating the 3-state kagome-lattice Potts antiferromagnet at zero temperature. We prove that this algorithm is not ergodic for symmetric subsets of the kagome lattice with fully periodic boundary conditions: given an initial configuration, not all configurations are accessible via Monte Carlo steps. The same conclusion holds for single-site dynamics.Comment: Latex2e. 22 pages. Contains 11 figures using pstricks package. Uses iopart.sty. Final version accepted in journa

"The Ising model on spherical lattices: dimer versus Monte Carlo approach"

We study, using dimer and Monte Carlo approaches, the critical properties and finite size effects of the Ising model on honeycomb lattices folded on the tetrahedron. We show that the main critical exponents are not affected by the presence of conical singularities. The finite size scaling of the position of the maxima of the specific heat does not match, however, with the scaling of the correlation length, and the thermodynamic limit is attained faster on the spherical surface than in corresponding lattices on the torus.Comment: 25 pages + 6 figures not included. Latex file. FTUAM 93-2

Dynamic Critical Behavior of a Swendsen-Wang-Type Algorithm for the Ashkin-Teller Model

We study the dynamic critical behavior of a Swendsen-Wang-type algorithm for the Ashkin--Teller model. We find that the Li--Sokal bound on the autocorrelation time ($\tau_{{\rm int},{\cal E}} \ge {\rm const} \times C_H$) holds along the self-dual curve of the symmetric Ashkin--Teller model, and is almost but not quite sharp. The ratio $\tau_{{\rm int},{\cal E}} / C_H$ appears to tend to infinity either as a logarithm or as a small power ($0.05 \leq p \leq 0.12$). In an appendix we discuss the problem of extracting estimates of the exponential autocorrelation time.Comment: 59 pages including 3 figures, uuencoded g-compressed ps file. Postscript size = 799740 byte

Absence of Phase Transition for Antiferromagnetic Potts Models via the Dobrushin Uniqueness Theorem

We prove that the $q$-state Potts antiferromagnet on a lattice of maximum coordination number $r$ exhibits exponential decay of correlations uniformly at all temperatures (including zero temperature) whenever $q > 2r$. We also prove slightly better bounds for several two-dimensional lattices: square lattice (exponential decay for $q \ge 7$), triangular lattice ($q \ge 11$), hexagonal lattice ($q \ge 4$), and Kagom\'e lattice ($q \ge 6$). The proofs are based on the Dobrushin uniqueness theorem.Comment: 32 pages including 3 figures. Self-unpacking file containing the tex file, the needed macros (epsf.sty, indent.sty, subeqnarray.sty, and eqsection.sty) and the 3 ps file

Transfer Matrices and Partition-Function Zeros for Antiferromagnetic Potts Models. V. Further Results for the Square-Lattice Chromatic Polynomial

We derive some new structural results for the transfer matrix of square-lattice Potts models with free and cylindrical boundary conditions. In particular, we obtain explicit closed-form expressions for the dominant (at large |q|) diagonal entry in the transfer matrix, for arbitrary widths m, as the solution of a special one-dimensional polymer model. We also obtain the large-q expansion of the bulk and surface (resp. corner) free energies for the zero-temperature antiferromagnet (= chromatic polynomial) through order q^{-47} (resp. q^{-46}). Finally, we compute chromatic roots for strips of widths 9 <= m <= 12 with free boundary conditions and locate roughly the limiting curves.Comment: 111 pages (LaTeX2e). Includes tex file, three sty files, and 19 Postscript figures. Also included are Mathematica files data_CYL.m and data_FREE.m. Many changes from version 1: new material on series expansions and their analysis, and several proofs of previously conjectured results. Final version to be published in J. Stat. Phy

Spanning forests and the q-state Potts model in the limit q \to 0

We study the q-state Potts model with nearest-neighbor coupling v=e^{\beta J}-1 in the limit q,v \to 0 with the ratio w = v/q held fixed. Combinatorially, this limit gives rise to the generating polynomial of spanning forests; physically, it provides information about the Potts-model phase diagram in the neighborhood of (q,v) = (0,0). We have studied this model on the square and triangular lattices, using a transfer-matrix approach at both real and complex values of w. For both lattices, we have computed the symbolic transfer matrices for cylindrical strips of widths 2 \le L \le 10, as well as the limiting curves of partition-function zeros in the complex w-plane. For real w, we find two distinct phases separated by a transition point w=w_0, where w_0 = -1/4 (resp. w_0 = -0.1753 \pm 0.0002) for the square (resp. triangular) lattice. For w > w_0 we find a non-critical disordered phase, while for w < w_0 our results are compatible with a massless Berker-Kadanoff phase with conformal charge c = -2 and leading thermal scaling dimension x_{T,1} = 2 (marginal operator). At w = w_0 we find a "first-order critical point": the first derivative of the free energy is discontinuous at w_0, while the correlation length diverges as w \downarrow w_0 (and is infinite at w = w_0). The critical behavior at w = w_0 seems to be the same for both lattices and it differs from that of the Berker-Kadanoff phase: our results suggest that the conformal charge is c = -1, the leading thermal scaling dimension is x_{T,1} = 0, and the critical exponents are \nu = 1/d = 1/2 and \alpha = 1.Comment: 131 pages (LaTeX2e). Includes tex file, three sty files, and 65 Postscript figures. Also included are Mathematica files forests_sq_2-9P.m and forests_tri_2-9P.m. Final journal versio

Transfer matrices and partition-function zeros for antiferromagnetic Potts models. VI. Square lattice with special boundary conditions

We study, using transfer-matrix methods, the partition-function zeros of the square-lattice q-state Potts antiferromagnet at zero temperature (= square-lattice chromatic polynomial) for the special boundary conditions that are obtained from an m x n grid with free boundary conditions by adjoining one new vertex adjacent to all the sites in the leftmost column and a second new vertex adjacent to all the sites in the rightmost column. We provide numerical evidence that the partition-function zeros are becoming dense everywhere in the complex q-plane outside the limiting curve B_\infty(sq) for this model with ordinary (e.g. free or cylindrical) boundary conditions. Despite this, the infinite-volume free energy is perfectly analytic in this region.Comment: 114 pages (LaTeX2e). Includes tex file, three sty files, and 23 Postscript figures. Also included are Mathematica files data_Eq.m, data_Neq.m,and data_Diff.m. Many changes from version 1, including several proofs of previously conjectured results. Final version to be published in J. Stat. Phy